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1. Introduction 

In the revolutionary braneworld viewpoint, our universe is a 3-brane embedded in an 
extra dimensional bulk. Standard matter and all interactions are confined on the brane; 
only graviton and possibly non-standard matter are free to probe the full bulk pp. 
Based on the braneworld viewpoint, our universe may contain many more dimensions 
than those we experience with our senses. The most compelling reasons to believe in 
extra dimensions are that they permit new connections between physical properties 
of the observed universe and suggest the possibility for explaining some of its more 
mysterious features. Extra dimensions can have novel implications for the world we see, 
and they can explain phenomena that seem to be mysterious when viewed from the 
perspective of a three-dimensional observer. Even if one is doubtful about string theory 
due to, for instance, its huge number of landscapes, recent researches have provided 
perhaps the most compelling argument in the favor of extra dimensions: a universe 
with extra dimensions might contain clues to physics puzzles that have no convincing 
solutions without them. This reason alone makes extra dimensional theories worthy of 
investigation. In this streamline, the braneworld models that are inspired by ideas from 
string theory provide a rich and interesting phenomenology, where higher- dimensional 
gravity effects in the early and late universe can be explored, and predictions can be 
made in comparison with high-precision cosmological data. Even for the simplest models 
of Randall- Sundrum (RS) [2] and Dvali-Gabadadze-Porrati (DGP) [3], braneworld 
cosmology brings new implications on the inflation and structure formation. Also it 
brings new ideas for dark energy and opens up exciting prospects for subjecting M- 
theory ideas to the increasingly stringent tests provided by high-precision astronomical 
observations. At the same time, braneworld models provide a rich playground for 
probing the geometry and dynamics of the gravitational field and its interaction with 
matter. In these respects, the DGP braneworld model is a scenario that gravity is 
altered at immense distances by the excruciatingly slow leakage of gravity off our 3- 
brane universe. In this braneworld scenario, the bulk is considered as empty except for 
a cosmological constant and the matter fields on the brane are considered as responsible 
for the evolution on the brane [3]. The self-accelerating DGP branch explains late-time 
speed-up by itself without recourse to dark energy or other mysterious components p|]. 
Even the normal DGP branch has the potential to realize an effective phantom phase 
via dynamical screening of the brane cosmological constant [5]. 

Scalar fields play a crucial role in modern cosmology, both in models of the early 
universe and late-time acceleration. Scalar fields provide also a simple dynamical model 
for matter fields in a braneworld and dark energy models. In the early universe, inflaton 
as a scalar field provides the required basis of the some well-established inflation models. 
Also at late time, dark energy models based on dynamical scalar fields have been studied 
extensively in recent years [6]. In braneworld models, the existence of a scalar field on 
the brane provides a variety of possibilities that brings the corresponding theory to 
explain some novel properties. In fact, a particular form of the bulk or brane matter is 
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a scalar field. In the context of braneworld induced gravity, it is natural to consider a 
non-minimal coupling of the scalar field and induced Ricci curvature on the brane. The 
resulting theory can be thought of as a generalization of the BransDicke type scalar- 
tensor gravity in a braneworld context [7J . As has been pointed in [8] , the introduction of 
the non-minimal coupling (NMC) is not just a matter of taste: the NMC is forced upon 
us in many situations of physical and cosmological interest. For instance, NMC arises at 
the quantum level when quantum corrections to the scalar field theory are considered. 
Even if for the classical, unperturbed theory this NMC vanishes, it is necessary for the 
renormalizability of the scalar field theory in curved space. In most theories used to 
describe inflationary scenarios, it turns out that a non-vanishing value of the coupling 
constant is inevitable. In general relativity, and in all other metric theories of gravity 
in which the scalar field is not part of the gravitational sector, the coupling constant 
necessarily assumes the value of | jH]- Therefore, it is natural to incorporate an explicit 
NMC between the scalar field and the Ricci scalar in the inflationary paradigm and 
also in scalar fields models of dark energy. In particular the effect of this NMC in a 
DGP-inspired braneworld cosmology has been studied by some authors (see [7J and also 

S)- 

There are several studies focusing on braneworld models with brane/bulk scalar 
fields. Some of these studies concentrate on the bulk scalar fields minimally or 
nonminimally coupled to the bulk Ricci scalar [TO]. The scalar field minimally or non- 
minimally coupled to gravity on the brane are studied by some authors [TTJ [TJl [TOJ [H] . 
In [2], the authors are studied the self-accelerating solutions in a DGP brane with a 
scalar field trapped on it within a dynamical system perspective. They have shown that 
the dynamical screening of the scalar field self-interaction potential occurring within 
the Minkowski cosmological phase of the DGP model mimics 4D phantom behavior 
and is an attractor solution for a constant self-interaction potential. But, this is not the 
case necessarily for an exponential potential. For exponential potential, they have shown 
that gravitational screening is not even a critical point of the corresponding autonomous 
system of ordinary differential equations. Along with this pioneer work, we consider a 
scalar field trapped on the DGP brane and we suppose this scalar field is non-minimally 
coupled to induced gravity on the brane. We study cosmological dynamics on the 
normal branch of the scenario within a phase space approach with both quintessence 
and phantom fields on the brane. We provide a phase space analysis of each model 
through a detailed study of the fixed points, their stability and cosmological viability 
of the solutions. We also study the classical stability of the solutions in each case in 
the w v — w' phase-plane. Our study, in comparison with existing literature in this 
field, provides a complete framework and contains several aspects of the problem not 
been considered yet. Since the self-accelerating DGP branch has ghost instabilities, our 
study here is restricted just to the normal DGP branch of the models. While the normal 
branch of a pure DGP setup has not the potential to explain the late-time cosmic speed- 
up and crossing of the phantom divide, we show that with a scalar field on the brane 
there are several new possibilities in the favor of these observationally supported issues. 
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We note that our motivation to study this model is as follows: as we have indicated 
above, observations support (at least mildly) that the equation of state parameter of 
dark energy has crossed the cosmological constant line (w = —1) in recent past (at 
redshift z ~ 0.25). It is impossible to realize this feature with a quintessence or phantom 
field minimally coupled to gravity in standard 4-dimensional theory [6j. Although the 
original DGP model was proposed to realize accelerated expansion of the universe in 
a braneworld setup, the self-accelerating branch of the DGP cosmological solutions 
has ghost instability. It is impossible also to cross the phantom divide line without 
a scalar field in the self-accelerating DGP branch [II] . The normal DGP solution 
has no ghost instability, but it cannot explain accelerated expansion and crossing of 
the phantom divide. It has been shown that localizing a scalar field on the normal 
DGP setup realizes these features [H]. It is possible also to incorporate extra degrees 
of freedom on the braneworld setup to have more successful models (see for instance 
[TJ. These extra degrees provide new facilities and richer cosmological history on the 
brane, a part of which is related to the wider parameter space. On the other hand, 
considering just a cosmological constant on the brane, although explains accelerated 
expansion through dynamical screening of the brane cosmological constant, it has not 
the potential to explain crossing of the phantom divide [5] . In this paper we have shown 
that a scalar field, minimally or nonminimally coupled to induced gravity on the brane 
has the potential to fill these gaps. We stress that all of the accelerated phases obtained 
in this study belong to the normal DGP branch of the model. Our study, based on the 
phase space analysis, provides the most complete treatment of the issue in the field. We 
have provided a complete analysis of the problem focusing on all possible details. 

2. Cosmological dynamics of a minimally coupled scalar field on the DGP 
brane 

DGP braneworld scenario has attracted a lot of attention through these years. Although 
this scenario has very interesting phenomenological aspects, it suffers from shortcomings 
such as ghost instabilities in its self-accelerating branch of the solutions. Fortunately 
the normal, non-self-accelerating branch of the DGP cosmological solutions has no ghost 
instabilities. For this reason we consider a scalar field on the DGP brane and we analyze 
cosmological dynamics of the normal DGP solutions in this setup. To begin and in 
order to explain the frame of our analysis, we start with the case of a minimally coupled 
quintessence field on the normal DGP branch. We note that this problem has been 
considered previously in [1]} [12], [13]. Especially, in [H] the authors have presented a 
detailed study of the problem with just a minimally coupled scalar field on the brane. 
We firstly study quintessence field with more details and more enlightening analysis 
than [13] both in calculations and corresponding analysis of the phase space points. 
Then we extend our study to the minimally coupled phantom fields, nonminimally 
coupled quintessence fields and finally non-minimally coupled phantom fields on the 
brane. In each step, we provide a detailed analysis of the model in phase space and 
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within a dynamical system approach. We study the late-time cosmological viability of 
the solutions and their stability in each step. We also investigate the classical stability 
of the solutions in w v — w' phase-plane. 

2.1. Minimally coupled quintessence field on the brane 

A minimally coupled quintessence field on the normal DGP within a phase space 
approach first has been considered in [TT]. There, the authors have considered the 
cosmological evolution of a QDGP model in a phantom-like prescription. They have 
defined an effective cosmological constant that increases by time evolution of the Hubble 
parameter and therefore realizes a phantom-like behavior through dynamical screening 
of the brane cosmological constant. Here we adopt another strategy: we focus on the 
phase space coordinates instead of the brane cosmological constant similar to strategy 
adopted in [14J . Although in this step the results are the same, but our analysis in this 
subsection introduces the notation, conventions and general framework of our procedure. 
In addition, some novel ingredients such as — w' phase-plane stability analysis, more 
detailed calculations and more enlightening plots are presented. 

The action of an induced gravity braneworld model can be written as follows (we 
use the sign convention of [6J) 

~ M i f j5 v r~„v & 



S 



■ J d 5 X^/=g~R 5 - y / d 4 xV^h~R<i + J d^xV^hC m + S GH , (1) 

where g a b is the metric of the bulk manifold and h^ is the induced metric on the 
brane. R 5 and R4 denote the 5 and 4 dimensional Ricci scalars respectively, and C m is 
the matter Lagrangian confined on the brane. Sqh is the Gibbons-Hawking boundary 
action which is required in order to apply the boundary conditions properly. In this 
induced gravity braneworld setup, the ratio of the two scales, the 4-dimensional Planck 
mass /1 and its 5-dimensional counterpart M 5 , defines the DGP crossover scale as 



2M| 



which determines the behavior of gravity in different distance scales on the brane. 
Adopting a FRW line element, the cosmology of the model is based on the following 
Friedmann equation |4l |TT] 

^ + f = (\Z^ + 4^ + ^) ' (2) 
where p is the energy density of the total cosmic fluid on the brane and consists of the 
energy densities of the scalar field and ordinary matter on the brane. We consider the 
flat geometry in which (K = 0), so the above equation reduces to 



» 2 + - = ^, (3) 
r c 3/i 2 



for the normal DGP branch of the scenario. When the Hubble length if 1 is much 
smaller than r c , which stands for the early time, the term — can be ignored relative 
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to the first term on the left hand side of ([3D- This term becomes important on scales 
comparable to the crossover scale when H' 1 is larger than r c , which corresponds to 
the late-time of the universe evolution. We emphasize that in which follows we focus 
just on the normal, non-self-accelerating DGP branch of the solutions. This is because 
the normal branch has no ghost instabilities. Nevertheless, the pure normal branch 
cannot explain the late-time acceleration without additional components on the brane. 
However, in the presence of a scalar field on the brane, it is possible to realize the late- 
time cosmic speed-up even in the normal DGP branch [TT1 [Ej. This is the reason why 
we considered an extension of the DGP setup with a scalar field on the brane. 

By differentiation of fl3]) with respect to the cosmic time, we obtain an equation for 
evolution of the Hubble parameter 

- {p+P) ( 1 1 ^ 
2/i 2 V 1r c H< 

where we have used the continuity equation on the brane as p + 3H(p + p) = 0. We 
note that the negativity of H ensures the phantom-like behavior on the brane The 
dynamics of the scalar field localized on the brane is given by the following Klein-Gordon 
equation 

dV 

V + 3Htp = -—. 5 
dip 

The energy density and pressure of the total matter localized on the brane are given by 
P = P v + p m = \<P 2 + V(<p) +p m , (6) 

P = P<p +Pm = l^P 2 ~ V(<p) + WmPm (7) 

respectively. To translate our equations of the cosmological dynamics in the language of 
the autonomous dynamical system, we introduce the following dimensionless quantities 



if vv y/p. 

x l ~ ~~r^ — 77 > x 2 ~ ~~r^ — 77 > x 3 



y/QpH ' y/3pH ' y/3pH ' /q\ 

1 r - VV " 

where a prime marks differentiation with respect to the scalar field, /=■£-. Through 
this paper we consider the exponential potential of the scalar field as V(ip) = V e~ Xip ^. 
This potential corresponds to a constant A and gives T = 1. The case of a constant 
potential (as has been studied separately in [2]), is a special case of this potential. We 
note that the above dimensionless quantities have explicit physical origin: x\ is related 
to the kinetic energy of the field, x\ is related to the potential energy of the scalar field, 
x\ is related to the ordinary matter density on the brane, x\ reflects the DGP character 
of this setup, and A and T are actually the slow-roll parameters of the model. Our main 
equations are the Friedmann equation that appears as a constraint, the Klein-Gordon 
equation and the continuity equation on the brane. These equations with the above 
dimensionless quantities provide the basis of our dynamical system analysis. 
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By using the above dimensionless quantities, we rewrite (EJ), flSJ) and ([7]) in the 
following form 

(p = -3fiH 2 (V6x 1 - Xx\) , (9) 

p = 3fi 2 H 2 ( y x 2 1 +x 2 2 + x 2 3 ), (10) 

p = 3p 2 H 2 (x\ — x\ + w m x\) . (11) 
The effective equation of state parameter in this case is given by 

9 9 9 

rp ^ ry* ^ _l_ /in rjr> & 

Weff = x 2 + x 2 + x, 2 ■ (12) 
We can obtain the exact cosmological solutions at the critical points by using the 
following relation 

■ _ 3 [2x 2 + (l + w m )x 2 ] 2 

2 1 + xl 1 1 

In each fixed point, we can rewrite this relation in the following form 

H = --H 2 , (14) 
a 

where by definition 



3 



3(2x 1 2 + (l + w m )x 

An integration of ( TT4"1) with respect to the cosmic time gives 

a(t) =a (t-t ) a (16) 

which corresponds to an accelerating phase if a > 1. Note that equation (16) is valid 
only in a small neighborhood around the critical point where a can be considered to be 
nearly constant. Now the Friedmann constraint equation in terms of the dimensionless 
quantities becomes 

xj + xj + xj- 2x 2 4 = l. (17) 

So, the allowable region of the phase space is actually outside of a unit sphere defined 
as x\ + x\ + x\ = 1. Now the autonomous dynamical equations are given as follows 

dx l o , V^, 2 /2X! 2 + (1 + W m )x£\ 

— — = -3x 1 + -—Xxi + 3xi t. 5 \ , (18) 

dN 2 2 \ 1 + x? + x 2 2 + x£ J K ' 



dx 2 VQ (2x x 2 + (1 + w m )xf « 

T77 = — —\xix 2 + 3x 2 5 5 |- , (19) 

diV 2 V 1 + x? + x 2 2 + x 2 1 



c^3 3 / 2x x 2 + (1 + w m )x 3 2 ^ 

rfiV = "2 (1 + " m)X3 + 3X3 ( l + x 1 2 + x 2 2 + x 3 2 J ' (20) 
where by definition, N = lna(t). The eigenvalues of the Jacobian matrix are as follows 
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Table 1. Location and dynamical character of the fixed points. 



name 




x 2c 


^3c 


Existence 


stability 




If 




«(() 


(la), (lb) 








±1 


VA, 7 


saddle point 





undefined 


7 - 1 


ao(t-t ) 2/3T 


(2a) 


1 








VA, 7 


saddle point for A > y/6 
unstable node for A < V*6 


1 


2 


1 


ooC - *o) 1/3 


(2b) 


-1 








VA, 7 


unstable node for A > — y/6 
saddle point for A < — \/6 


1 


2 


1 


ao(t - t ) 1/S 


(3a), (3b) 


A 

y/S 


±(1"^) 1/2 





A 2 < 6 


saddle point 


1 


3 


^-1 
3 1 


ao(* - «o) 2/A ^ 


(4a), (4b) 
(4c), (4d) 


fSl 

y 2 x 




±(1-^)V2 


A 2 > 37 


stable node for 37 < A 2 < gzrz '2 
stable spiral for A 2 > 2 


A* 


7 


7 - 1 


a (t - t ) 2/3T 



points la , lb : 



3 3 3. 

"i = 2~7, «2 = "3 = -2( 2- 7) • 



point 2a : 



ai = 3, a 2 = — |^A + 3, a 3 = ^(2-7) 



point 2b : 
points 3a , 3b : 



at = 3, a 2 = ^2~A + 3, a 3 =~(2-7). 



«i = y, «2 = y-3, a 3 = i(A 2 - 37) . 



points 4a , 4b , 4c , 4d : 
3 



ai = ^7, 



«2,3 




S7(A 2 -37)' 



A 2 (2- 7 ) 



where 7 = 1 + w m is the barotropic index which depends on the type of ordinary matter 
on the brane. In what follows, we consider < 7 < 2. 

Table [U summarizes the complete information about existence, stability and 
cosmological characteristics of these phase points. In this table, = y- where 
p c = 3fi 2 H 2 and 7^ = 1 + w v . We note that for point 1 with 7 = 0, there is a 
stable center. All other mentioned points have center too, depending on the values of 7 
and A. But these centers are not necessarily stable. For instance, for point {(3a), (3b)}, 
the center is unstable. 

For points {(la), (lo)}, there is no contribution of the scalar field and the universe 
is dominated by other matter fields. According to the eigenvalues and since < 7 < 2, 
these points behave like saddle points in the phase space. 

Points {(2a), (26)} are solutions dominated by the kinetic energy of the scalar field. 
The contribution of scalar field potential and the energy densities of other matter fields 
are irrelevant for these phases. For these points, we obtain 7^ = 2, which is referred to 
as a stiff matter equation of state. These points have no late time acceleration and the 
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Table 2. Location and critical point for 7^ = (that is, w„ = — 1). 



name 


%2c 


%3c 




stability 


W eff 


a(t) 


(la),(lb) 





±1 





unstable nodes 


7-1 


(t - tof^ 


(2a),(2b) 


±(l + 2x 4 2 ) 1 /2 





£4 


stable attractor 


-1 


e A(t-to) (A = constant) 



stability of them depends on the values of A, so that for A < (for point (2a)) and 
A > — \/6 (for point (26)), we obtain unstable nodes. Otherwise they are saddle points. 

For critical points {(3a), (36)}, the energy density of universe is dominated by the 
scalar field's kinetic and potential energies. For A 2 < 2 we have accelerated phase 
of expansion, but this phase is not stable. In these cases, for A = 0, the universe is 
dominated by a cosmological constant. 

The last line of table Q] contains four critical points {(4a), (46), (4c), (4a 1 )}, so that 
depending on the values of 7 and A, we have stable spirals or stable nodes. Here 7^ 
is equal to the barotropic index of matter, 7. This is a reflection of the fact that the 
exponential potential used in this framework can give rise to an accelerated expansion 
and possesses cosmological scaling solutions in which the field energy density p v is 
proportional to the fluid energy density, p m . 

The ingredients of table [1] can be explained with more geometrical details through 
the phase space trajectories. Figure [U shows the two dimensional (x\ — X2) phase plane 
for A = +1 . In this figure, points A and B both are unstable nodes. 

We note that the above arguments were based on the assumption that x 4 = for 
all critical points. Since X4 is related directly to the braneworld nature of the solutions, 
our analysis up to this point was effectively 4- dimensional. Now we consider the case 
that X4 7^ 0. Table [2] gives the results of the corresponding phase space analysis. In 
this case, there is a cosmological constant dominated accelerating phase which is an 
stable attractor corresponding to critical lines (2a), (26) of table [2] (curves C\ and C2 of 
figure |5J) • 

As we have stated previously, the pure normal DGP solution has not the potential to 
explain the late-time acceleration. In our case with a quintessence field on the brane, as 
we have shown, it is possible to realize the late-time acceleration in this setup. We note 
also that as table [Q shows, for fixed points {(la), (lb)}, {(3a), (3b)} and also {(4a), (4b)} 
it is possible in principle to realize a late-time acceleration which depends on parameters 
7 and A and the stability of which needs to be investigated in each case. For fixed points 
{(la), (lb)} (first row of tabled]), the accelerating phase (with q < where q = —ft) is 
possible if 7 < 2/3. But this accelerating phase is a saddle point and obviously is not the 
late-time cosmic accelerating phase. For fixed point {(3a), (3b)}, the accelerating phase 
is possible if A 2 < 2. This gives also a saddle point which is not corresponding to the late- 
time stable, accelerating phase of cosmic expansion. For fixed point {(4a), (4b), (4c), 
(4d)}, there is an accelerating, stable phase if 7 < 2/3. However, these stable points are 
either a node or a spiral. Albeit these are not corresponding to a late-time accelerating, 
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Figure 1. The phase plane for A = +1 . Points A and B are unstable nods. 




Figure 2. The phase plane for w m = 0. Point A is an unstable node(which 
reflects the first line of tabled]). There are two curves {C\ and C2), which are 
corresponding to critical lines ((2a), (26)) of table [2j 



stable, de Sitter attractor. So, none of the accelerating phases corresponding to tabled] 
are the de Sitter phase. Nevertheless, as we have stated previously, point (2) (actually 
a critical line) of table [2] gives a de Sitter attractor. Therefore, with a quintessence field 
in the normal DGP setup, it is possible to realize the late-time acceleration in contrast 
to pure DGP case. 

Now we analyze the accelerating phase of the model through the evolution of the 
deceleration parameter. The equation for the deceleration parameter is q = — |f = 
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|(1 + 3w e ff). In our case, it can be written by using the phenomenological parameters 



Vt m and Vtp as 



1 (1 + Sw m ) Q m (1 + z 



3(l+«i m ) 



+ (i + 3w v )n v (i + z 



.3(l+w 



<p) 



2 a m (i + 2 ) 3(1+ " m) + ^(i + ^) 3(1+ ^ ) 

Figure [3] shows the behavior of g versus the redshift. There is a transition to the 
accelerating phase at z — 0.68. 



-O.l 




Figure 3. The deceleration parameter versus the redshift for fi r 



0.28, 



= 0.8, w r , 
at z = 0.68. 



and w v = —2/3. Transition to the accelerating phase occurs 



We define the dimensionless density parameters as follows 

Po PtpO ~ 1 



3ff n 2 ' 



3ff n 2 ' 



Then the Friedmann equation can be rewritten as 



if (z) = ff ( y/n rc + Q m (1 + z f {1+w ™ ] + fi„ (1 + z) 3(1+ ^ 



- y^T) (22) 



Figure H] shows the evolution of the Hubble parameter versus the redshift and equation 
of state parameter of a minimally coupled quintessence scalar field with Q m = 0.28, 
Q v = 0.8. We can understand from this figure that the Hubble parameter of the model 
decreases as the redshift decreases. This is a trace of essentially possible realization of 
an effective phantom-like behavior on the brane (see [Tlj for more details). Figure [5] is 
a 2-dimensional plot of if versus the redshift for a quintessence field with = —0.8. 

Another important issue is the possibility of crossing of the cosmological constant 
line (the so-called phantom-divide line) by the equation of state parameter. The effective 
equation of state parameter can be written as follows 

Wm n m (i + + Wv n v (i + z f l+w * ] 

w eff — n ,„ , ,3(1+%) ~ /- v an -4- ,«..-> • u "" 



,3(1+ «v) 
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Figure 4. The 3-dimensional plot of the Hubble parameter versus the redshift 
and equation of state parameter of the scalar field. 




Figure 5. This figure shows the plot of the Hubble parameter versus the 
redshift for a quintessence field. It is plotted for Q m = 0.28, = 0.8, w m = 
and Win = —0.8. 



As figure [6] shows, it is impossible to cross the cosmological constant equation of state 
parameter w v = —1 by a minimally coupled quintessence field in the normal DGP setup. 

Now we focus on the classical stability of the solutions in — w' phase-plane 
of the present model (see [15] for a similar analysis for other interesting cases). Here 
a prime marks the derivative of w v with respect to the logarithm of the scale factor, 
N = lna(t), so that 

= dw^ = dw^dp^ 
v dN dp^ dN' 1 ' 

We define the function c a so that c\ = 4^ or equivalently c\ = j^. Generally the 
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-2-10 1 2 



Figure 6. The effective equation of state parameter versus the redshift for 
U m = 0.28, 0^ = 0.8, w m = and w v = —2/3. There is no crossing of the 
cosmological constant line. 



sound speed expresses the phase velocity of the inhomogeneous perturbations of the 
scalar field. If we suppose the scalar field's energy-momentum to have a perfect fluid 
form, this function would be the adiabatic sound speed of this fluid. To avoid the future 
big rip singularity, we set c\ > 0. Since 



(25) 



and 



we obtain 



dvjy _ J_dp^ _ 
dp v p v dp v p£ 
1 f dp v 

= — b — w v 

p v K dp v 

^ = i = -3(l + -,K, (26) 
-3(1 + w v ) (c 2 a - . (27) 



Therefore, we obtain 

v/ v = -3(1 - w*) + A^/3(l + ^)^(l - w v ) . (28) 
Now the w v — w'^ phase plane is divided into the following two regions 

w v > -1 , w'p < 3w^(l + w v ) c\ > (region I) 

(29) 

w ip > —I ,w' tp > 3w v (l + w^) cl < (region II) 

that are shown in figure [71 The region I is the classical stability region of the theory. We 
note that there is no phantom phase in this case and therefore we have not encounter 
with four distinct regions of — w' plane as usually are discussed in literature. 
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Vf<p 

Figure 7. Bounds on w' as a function of in w v — w' phase-plane for 
Sl v = 0.8 and A = 0.1 

2.2. Minimally coupled phantom field on the normal DGP setup 

Astrophysical data indicate that w, the equation of state parameter of the cosmic fluid, 
lies in a very narrow strip close to —1. The case w = — 1 corresponds to the cosmological 
constant. For w less than —1, the phantom dark energy is observed, and for w more 
than —1 (but less than— ~) the dark energy is described by a quintessence field which 
has been studied in the previous subsection. Moreover, the analysis of the properties 
of dark energy from recent observational data mildly favor models of dark energy with 
w crossing the —1 line in the near past. So, the phantom phase equation of state with 
w < — 1 is mildly allowed by observations [161 El HH] • I n this case, the universe currently 
lives in its phantom phase which ends eventually at a future singularity (the Big Rip 
singularity [19J). There are also a lot of evidence all around of a dynamical equation of 
state, which has crossed the so called the phantom divide line w = — 1 recently, at the 
value of redshift parameter z ~ 0.25 [161 El HE]- For these reasons, now we pay our 
attention to a phantom field localized on the DGP brane and we study cosmological 
dynamics of the normal DGP branch in this case within a phase space analysis. 



The energy density and pressure of a phantom field are defined as 

Pv = -\<P 2 + V(<p) , (30) 

V* = ~ V(<p) , (31) 

respectively. The Klein-Gordon equation governing on the dynamics of the phantom 
field is given by 

dV , s 

(p + 3Hp = — . (32) 

dip 
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Table 3. Location and dynamical character of the fixed points. 



name 




x 2c 


^3c 


Existence 


stability 




If 


«>eft 


a(t) 


(la), (lb) 








±1 


VA, 7 


saddle point 





undefined 


7 ~ 1 


«o(< - to) 2737 


(2a), (2b) 


A 


±(1 + ^) 1/2 





VA, 7 


stable node 


i 


-a" 

3 


^-1 

3 


oo(< - to) _2/Ai 


(3a), (3b) 
(3c), (3d) 


y 2 x 


±( 3 2 |^ )1/2 


±(l + l^) 1/2 


7 < , A 2 > -37 


saddle point 


3 T 


7 


7 - 1 


2 

O0(< - *0) 3 7 



With phase space coordinates as defined in (jSJ), The Friedmann equation of the 
model in phase space now is written as 

- x\ + x 2 + x 2 - 2x 2 A = 1 > (33) 

and the autonomous dynamical equations are 



etei \/6. o ( -2x? + (I + w)xn i 

_ = -3xi - -7^ Ax 2 , + 3xi — 9 9 » 34 

diV 2 2 1 1 - xi 2 + x 2 2 + x 3 2 1 



dx 2 y/6 + (1 +W)X 3 2 ^ 

377 = — ^Axix 2 + 3x 2 2— 2— 2" , (35) 

dN 2 ^ 1 — x* + x{ + x£ ) 

dx^ _3 / -2x 1 2 + (l + w)x 3 2 \ 

diV 2 X3 l N l-x 1 2 + x 2 2 + X3 2 J • 1 j 

In this case there are eight critical points that are shown in table [3j The eigenvalues 
of these points are 

• points (la) , (lb): 

• points (2a) , (2b): 

11 A 2 
"i = -^(A 2 + 6), «2 = ~2(A 2 + 37), « 3 = -y 

• points (3a) , (3b) , (3c) , (3d): 



3 3 / 8t(P + 3t)\ 

ai = 2 7 ' «2,3 = -4(2-7)^l±^l- A2(2 _ 7) j- 

In table [3] we summarized also the results of the phase space analysis of the model in 
addition to cosmological characters of each critical point. For critical points {(la), (lo)}, 
there is no contribution of the phantom scalar field and the universe is dominated by 
matter fields other than the phantom scalar field. These points behave like saddle 
points in the phase space and for 7 < 2/3, these points give an accelerating phase. 
It is possible to have scaling solutions in this case too. For points {(2a), (26)}, there 
is no contribution of ordinary matter fields and the energy density of the universe is 
dominated by the phantom scalar field's kinetic and potential energies. In these cases, 
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there is no possibility of accelerated expansion on the brane. The corresponding points 
in phase plane are stable nodes. The last line of table [3] consists of four critical points 
{(3a), (36), (3c), (3c?)}. These are just saddle points in phase plane. There are scaling 
solutions in these critical points. For these fixed points accelerated expansion is possi- 
ble if 7 < 2/3. However, this accelerated expansion phase is not a de Sitter stable phase. 

We note that the analysis presented in the previous paragraph was based on the 
condition £4 = which gives essentially an effective 4-dimensional picture of the model. 
Now we consider the case that £4 7^ 0. The phase space analysis of the model with 
minimally coupled phantom field gives the same results as are presented in table [2] for a 
minimally coupled quintessence field. Similar to minimal quintessence scalar field case, 
in this case there is a cosmological constant dominated accelerating phase which is an 
stable attractor corresponding to critical lines (2a), (26) of table [2] (curves C\ and C 2 of 
figure [2]). So, with a phantom field minimally coupled to induced gravity in the normal 
DGP setup it is possible to have an attractor, de Sitter solution realizing the late-time 
accelerated expansion (see table [2]). Figure [8] shows a plot of the X\ — x 2 phase plane of 
the model. Point A is corresponding to points (la) and (lb) of table [3] Points B and C 
are corresponding to points (2a) and (2b) of table [3] which are stable nodes. The £3 — x 2 
(with £ 4 7^ 0) phase plane of the model is shown in figure [9j We note that this figure 
is actually the same as figure [2] but now plotted in £3 — x 2 plane rather than x 2 — £4 
plane. Here we encounter a line of stability points as is shown in figure [9] Points A and 
B are corresponding to critical points (la) and (16) of table EJ 

The form of the deceleration parameter and effective equation of state parameter 
for phantom field are the same as quintessence field presented as equations (21) and 
(23). Figure [10] shows the behavior of the deceleration parameter q(z). The universe 
enters the accelerated phase at z ~ 0.68. Also Figure [TT1 shows how the normal branch 
Hubble parameter evolves on the brane with a phantom field. Also Figure [12] shows 
the behavior of w e ff(z). There is a crossing of the phantom divide line in this setup. 
Therefore, with a minimally coupled phantom field on the normal DGP setup, it is 
possible to cross the phantom divide line by the effective equation of state parameter of 
the dark energy. Note that as we have shown previously, this crossing was impossible 
with a quintessence field on the brane. 

To investigate the classical stability of the solutions with phantom field in the 
normal DGP setup and within the w v — w'^ phase-plane approach, we adopt the same 
strategy as has been done for quintessence field. For this minimally coupled phantom 
field, we have 




(37) 



Now the w v — w' phase-plane is divided into the following two regions 

[ w v < -1 , w' < 3^(1 + w v ) c\ > (region I) 




(38) 



c^ < (region II) 
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Figure 8. The phase plane for A = +1 and w m = 0. Point A is a saddle 
point, whereas points B and C are stable nodes . 
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Figure 9. The X3 — £2 phase plane for w m = 0. Points A and -B are unstable 
nodes (which reflects the first row of table 2). There exists a critical line in this 
case. 



The Region I is the subspace of the classical stability of the solutions with a minimally 
coupled phantom field on the brane. 

3. Cosmological dynamics of a non-minimally coupled scalar field on the 
normal DGP setup 

Now we consider a non-minimally coupled scalar field on the DGP brane. As we have 
stated previously, in a realistic gravitational or cosmological scenario with scalar fields, 
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Figure 10. The deceleration parameter versus the redshift for Q m = 0.28, 
Vttp = 0.8, w m = and w v = —1.5. Transition to the accelerating phase occurs 
at z = 0.67. 
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Figure 11. Evolution of the Hubble parameter with redshift on the brane 
with a phantom field. Here we considered = —1.5. 

incorporation of the non-minimal coupling is inevitable. In fact, incorporation of a non- 
minimal coupling (NMC) between matter field and gravity is necessary from several 
compelling reasons. There are many theoretical evidences that suggest incorporation of 
an explicit non-minimal coupling of the scalar field and gravity in the action [201 |2T] . 
A nonzero non-minimal coupling arises from quantum corrections and it is required 
also for renormalizability of the corresponding field theory. Amazingly, it has been 
proven that the phantom divide line crossing of the dark energy described by a single, 



Stability analysis of the cosmological solutions with induced gravity and scalar field on the branel9 




t 1 1 1 1 1 1 1 r 

-2-10 1 2 



Figure 12. The effective equation of state parameter versus the redshift for 
f2 m = 0.28, $7^ = 0.8, w m = and w v = —1.5. Transition to the phantom 
phase occurs at z ~ 0.25. 
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Figure 13. Bounds on w' v as a function of w v in w v — w' v phase-plane for 
flip = 0.8 and A = 0.1. The subspace shown by Region I of the model parameter 
space is the classical stability subspace of the solutions. 

minimally coupled scalar field with a general Lagrangian is even unstable with respect 
to the cosmological perturbations realized on the trajectories of the zero measure [22J. 
This fact has motivated a lot of attempts to realize crossing of the phantom divide 
line by equation of state parameter of a scalar field non-minimally coupled to gravity 
as dark energy candidate in more complicated frameworks [23J. In which follows, we 
study cosmological dynamics with a non-minimally coupled quintessence/phantom field 
on the normal DGP setup within a phase space approach analysis. We study possible 
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realization of the late-time acceleration and crossing of the phantom divide line in this 
setup. We also investigate the classical stability of the solutions in separate regions of 
the w — w' phase-plane. 



3.1. Non-minimally coupled quintessence field on the normal DGP branch 

The equations governing on the cosmological dynamics on the normal DGP branch with 
a non-minimally coupled quintessence field follows 

H=-^i(^^r\ (39) 



2/i 2 V 2Hr c 

where p = p m + p v and p = p m + p v , and the energy density and pressure of the scalar 
field are defined as [7j 

Pv = l -p 2 + V(<p) + S£H<p<p + 3£#V , (40) 

and 

P v = \p 2 - V(<p) - 2^{ip(p + 2<pH<p + p 2 ) - & 2 (2H + 3H 2 ) , (41) 

respectively. As usual, p m is the energy density of ordinary matter fields other than the 
scalar field <p on the brane. The other dynamical equation is the following Klein-Gordon 
equation 

dV 

<p + 3Hcp + ZRcp=-—. (42) 

a<p 

Now we define the dimensionless variables as 



p vv 



(43) 



VEfiH' V3pH' V3pH' 

1 V? 
y/2r c H p 

By using these phase space variables, now the evolution equations ( 13"9|) and 
can be rewritten as 

• 3(1 - 20x x 2 + yoTxiXs + |7^ 3 2 + 3 v / !x 5 (v / 6x 1 + 4^x 5 - Xx 2 ) 

= l + xi-(l-60xi H ' (44) 

and 

-2^) (3(1 - 2Q Xl 2 + V^x lX5 + | 7 x 3 2 ) 



(p = -3pH 2 



1 + xi - (1 - 6£)x£ 

(1 + x? - 

+ 



1 + x 2 - x 2 )(VQx 1 + V£c 5 - Xx 2 ) 



1 + x 2 - (1 - 60x 5 2 



(45) 



respectively. Therefore, we obtain the following Friedmann constraint equation in the 
phase space of the model 

x 2 + x 2 2 + x 2 + 2^6^X1X5 + x 2 - 2x1 = 1 • ( 46 ) 
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To describe the dynamical system of the model, first we need to obtain the autonomous, 
phase space equations. We differentiate the phase space dimensionless variables with 
respect to iV = In a to find 

dx\ 



dN 

dx 2 
dN 
dx± 
dN 



—3xi + -^-Xx 2 2 



1 



and 



dx§ 
dN = 

where by definition, 



# = 



3(1 - 2£)x 2 + y^x lX5 + f 7 x 3 2 + 3^x 5 (VQx 1 + 4y/£x 5 - \x 2 2 ) 

l + xi-(l-QOxi 



(47) 
(48) 
(49) 

(50) 



The stability around the fixed points is related to the form of the eigenvalues in 
each critical point. The eigenvalues can be obtained by using the above autonomous 
equations, the results of which are as follows 



point (la), (lb): 



«1,2 = -1, "3,4 

point (2a), (2b): 



Q Q -1 

- + -7 ± -^36 - 367 + V - 192£ + 144£ 7 
2 4 4 



curve C: 



«i = --7 + 2, a 2 — —l, 0:3,4 = 2 



«i = --7> "2 = 0, 



- 48 f - 48 x 2 2 £ + 3 A 2 x 2 4 - 36 \ 2 x 2 4 £ 



1 1 

a3A = 2 (16 £ + 16 x 2 2 £ - \ 2 x 2 A + 12 A 2 x 2 4 _ 

± (2304 £ 2 + 4608 :r 2 2 £ 2 - 288 A 2 :r 2 4 £ + 1920 A 2 x 2 4 ^ 2 + 2304 x 2 4 £ 2 - 288 x 2 6 £ A 
+4992 x 2 Q i 2 \ 2 + 9 A 4 :r 2 8 - 72 A 4 :r 2 8 £ + 1872 A 4 :r 2 8 £ 2 - 24576 x 2 2 ^ - 12288 £ 3 
+384 A 4 x 2 6 £ - 3072 A 2 :r 2 2 £ 2 - 12288 x 2 ^ 3 - 12 X 6 x 2 10 - 9216 x 2 6 ^X 

-9216 A 2 x 2 4 ^ 3 + 144 \ 6 x 2 10 £ - 2304 A 4 :r 2 6 ^ * 
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Table 4. Location and dynamical character of the fixed points. 



name 


Xlc 


X2c 


X 3c 


X 5 c 


stability 


7 V 


W eff 


a(t) 


(la), (lb) 








±1 





saddle point 


undefined 


7-1 


'A 

ao(t - 1 ) 37 


(2a), (2b) 











±1 


saddle point 


4/3 


1/3 


ao{t - to)' 


C 





x 2 





4^ 


stable attractor 





-1 


e A(t-to) 



We note that when there is a zero eigenvalue for a critical point, it is necessary to 
use the center manifold theory in order to study the stability of that point in phase space 
of the model. In our case there is a zero eigenvalue for a critical line, and therefore there 
is no need to do the center manifold analysis. In other words, since we have a critical 
line here, the non-vanishing eigenvalues are enough in this case to treat the stability of 
the critical points (see for instance [21] )• Table H] summarizes the results of the stability 
analysis in the phase space of this model. Also this table contains types of possible 
cosmological dynamics in this setup. 

For critical points {(la), (16)}, there is no contribution of the scalar field and the 
universe is dominated by matter fields other than the quintessence scalar field. These 
two critical points behave like saddle points in the phase space. For 7 < 2/3, one can 
obtain an accelerating phase of expansion, but this phase is not stable. It is possible 
to have scaling solutions for these cases. As an interesting case, if 7 = 2/3, then we 
find w e ff = — 1/3 which shows domination of the curvature energy. The critical points 
{(2a), (26)} also behave like saddle points in the phase space and in these cases we have 
no late-time acceleration. For these points, the universe is radiation dominated. The 
last line of table HJ stands for a critical line (C) and in this case, there is a cosmological 
constant dominated accelerating phase. In this case, the potential energy of the brane 
scalar field plays the role of a cosmological constant on the brane. So, with a non- 
minimally coupled quintessence scalar field on the DGP brane, it is possible to realize 
a stable, de Sitter late-time accelerating phase even in the normal branch of the model. 
Figure [TH shows the phase plane of the model with £ = 1/6 (the conformal coupling). 
Point C is a stable attractor, whereas points A and B are saddle points. 

Figure [15] shows the three-dimensional, X\ — x 2 — £5 phase space of the model for 
the last row of table HI The stable, attractor points are located in a hyperbolic curve 
(curve C of table HJ) in x 2 — £5 perspective. In the three dimensional X\ — x 2 — x§ phase 
space, this is a hyperbolic hypersurface of stability points as shown in figure [TSJ All 
points of this hypersurface are stable attractors. We note that point C of figure [TJ] is 
corresponding to the mentioned hyperbolic on x\ — 25 plane. 

Figure [16] shows the possibility of phantom divide crossing by the equation of 
state parameter of the model. The universe transits into the phantom phase from a 
quintessence phase in a redshift that is observationally viable [z ~ 0.25). 
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Figure 14. The phase plane for £ 
points A and B are saddle points. 



1/6. Point C is a stable attractor, but 




From another perspective, the classical stability of the solutions in w v — w' v phase- 
plane gives some interesting results. To do this end, we calculate w v versus w'^ as 



< = 4^6f(l - 60(1 + w^lxuxs - 3(1 - 60^(1 + w^x 2 
-6 + 28£)x? + 2^(1 - 30xia; 2 2 - 8^(1 - A^) Xl x 5 + 2y^Xx 2 2 x 5 - 8£x 5 2 + 2 y /G£\ 2 x 1 x 2 2 x 5 

X X 2 + X 2 2 + X 5 2 + 2^/6£XlX 5 
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Figure 16. Crossing the phantom divide by the equation of state parameter 
for a(t) = a§e vt and <p = <poe~ at , where ao, </?o> v > and a > are constants. 
This crossing of the phantom divide occurs at z ~ 0.25 



-(x 2 +xi+x 2 +2^x 1 x 5 )b 3 (l-QO(l+w v ) 2 x 2 x 2 +3w^(l+w ip ) (l-b(l-6£)x B 



where by definition 



6 S (1 



2Hr, 



-i 



i + x: 



-i 



In this case, the w v — w'^ phase-plane is divided into the following four regions 
' w ip >-\, w' v> >3w ip {l + w tp ) ■■ 



w w <-\, w'>3wjl + w u 



w w > -1 , w' < SwJl + w u 



K w v <-l, w'< 3^(1 + 



c 2 < (region I) 

c 2 a > (region II) 

c 2 & > (region III) 

c 2 a < (region IV) 



(52) 



As we have explained in previous sections, the stability of the solutions requires c 2 > 0. 
So, the stability regions of the solutions in this case are the regions II and III. The 
region II corresponds to an effective phantom phase while region III is a quintessence 
phase. Figure [T71 shows these regions in w v — w' phase-plane. 

3.2. Non-minimally coupled phantom field on the normal DGP setup 

For completeness of our analysis, now we consider a non-minimally coupled phantom 
field on the normal DGP setup. The energy density and pressure of this non-minimally 
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Figure 17. Bounds on w' as a function of w v in w v — w' phase-plane for 



0.8, £ = 1/6 and A = 0.1 



coupled phantom field are given by 

1 



and 



p v = -^ 2 + V(<p) + 6£H<p<p + 3£#y , (53) 
2 V(ip) -2£(<p<p + 2<pH(p + <p 2 ) -& 2 {2H + 3H 2 ) . (54) 



2 
1 



Pv 



respectively. As previous cases, the equation governing on the the evolution of the Hub- 
ble parameter depends on the dimensionless variables and can be written as 



H 



3(-l - 2S)x 2 + ^Mxix 5 + hx 2 + 3^x 5 (VQx 1 - 4v^x 5 + A 



1 + xi - (1 + 60x 5 2 
The Klien- Gordon equation for a non- minimally coupled phantom field is 

dV 



H 2 . 



(p + 3H<p — ^Rip 



dip 



(55) 



(56) 



The Friedmann equation in terms of the phase space coordinates now can be written as 
the following constraint equation 

- x 2 + x 2 + xi + 2^6^x1X5 + x 2 - 2x1 = 1 • ( 57 ) 

The parameter space of the model as a dynamical system is described by the following 
autonomous system 

dx\ \/6 . 

dN = - 3Xi '^ 



dx2 
~dN 
dxi 
dN 



-3xi - — Xx 2 + 2^6^x5 + ( tJ6£x 5 - m)A 
V6. 



-\x1X2 — X2A . 



— X4A 



(58) 

(59) 
(60) 
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Table 5. Location and dynamical character of the fixed points. 



name 


Xlc 


X2c 


%3c 


%5c 


stability 


7 V 


w eff 


a(t) 


(la),(lb) 








±1 





saddle point 


undefined 


7-1 


'2 

a>o(t - 1 ) 37 


(2a),(2b) 











±1 


saddle point 


4/3 


1/3 




C 





x 2 





4^ 


saddle point 





-1 


e A(t-to) 



§ = V^*> ■ (61) 



where 



As previous sections, in order to study the stability of the critical points, we should 
obtain their eigenvalues, which are written as follows: 



point (la), (lb): 



point (2a), (2b): 



curve C: 



3 3 1 



"1,2 = 2-7, "3,4 = -- + -7 ± 4 V 36 - 36 7 + V + 192 £ - 144 ^7 



3 

«i = --7 + 2, a 2 = -1, «3,4 = 2 



3 

«i = --7, "2 = 0, 



"3,4 = r 



48 £ + 48 x 2 2 £ - 3 A 2 x 2 4 - 36 A 2 x 2 4 ^ 



2 (-16 £ - 16 x 2 2 £ + A 2 x 2 4 + 12 \ 2 x 2 4 

± (2304 £ 2 + 4608 x 2 2 £ 2 - 288 A 2 x 2 4 £ - 1920 A 2 x 2 4 £ 2 + 2304 x 2 4 £ 2 - 288 x 2 6 £ A 2 
-49 92 x 2 6 £ 2 A 2 + 9 A 4 x 2 8 + 72 A 4 x 2 8 £ + 1872 A 4 x 2 8 ^ 2 + 24576 x 2 2 £ 3 + 12288 £ 3 
-384 A 4 x 2 6 £ + 3072 A 2 x 2 2 ^ 2 + 12288 x 2 4 £ 3 + 12 A 6 x 2 10 - 9216 x 2 6 ^ 3 A 2 

-9216 A 2 x 2 4 £ 3 + 144 A 6 x 2 10 £ - 2304 A 4 x 2 6 £ 2 ^ 3 



Table O shows the results of the stability analysis in the phase space of a non- 
minimally coupled phantom field. In this case there are four critical points non of them 
result in a stable phase for the system. There is also a critical line, (C), which behaves 
as line of saddle points. 

Figure [TBI shows the x\ — £5 phase plane of the model with £ = 1/6. In this case 
points A, B and C behave as saddle points. 
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Figure 18. The x\ — X5 phase plane for a non-minimally coupled phantom 
field with £ = 1/6. Points A, B and C are saddle points. 
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Figure 19. Crossing of the phantom line by the equation of state parameter 
for a phantom field on the normal DGP setup. We considered a(t) = a§e vt 
and (p = ipoe~ at , where ao, , v and a are positive constants. Crossing of the 
phantom divide occurs at z ~ 0.24 



Now let us to study the possibility of phantom divide crossing by the equation 
of state parameter in this setup. As Figure [19] shows, in this case transition to the 
phantom phase occurs from quintessence phase at z ~ 0.25. Finally, investigation of 
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The Wp — w'p phase-plane of the model is plotted in figure 1201 This phase-plane can 
be divided into the following four regions 
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The stable regions of the phase-plane are those regions that the condition c 2 > is 
fulfilled. In this case, these stable regions are Region II and III as shown in figure | 




Figure 20. Bounds on w' v as a function of w v in w v — w' phase-plane for 
Clip = 0.8, £ = 1/6 and A = 1. The stability regions of the phase-plane are 
painted with cyan color. 



4. Summary and Conclusion 

The accelerated expansion of the universe supported by recent observational data could 
be associated with dark energy, whose theoretical nature and origin are still unknown 
for cosmologists. Cosmological constant or vacuum energy with an equation of state 



Stability analysis of the cosmological solutions with induced gravity and scalar field on the brane29 

parameter ou — — 1, is the most popular candidate for dark energy, but unfortunately it 
suffers from some serious problems such as huge fine-tuning and coincidence problems. 
Therefore, a number of models containing dynamical dark energy have been proposed 
as responsible mechanisms for late-time cosmic speed up. Some of these models are 
quintessence, k-essence, phantom scalar field, chaplygin gas models and so on. An- 
other alternative approach to explain the late-time cosmic speed up is modification of 
the geometric sector of the Einstein field equations leading to modified gravity theo- 
ries. In the spirit of modified gravity proposal, the Dvali-Gabadadze-Porrati (DGP) 
braneworld scenario explains the late-time accelerated expansion in its self-accelerating 
branch without need to introduce a dark energy component on the brane. However, 
some important features of dark energy such as possible crossing of the cosmological 
constant equation of state parameter are missing in the pure DGP model. In addition, 
the self-accelerating DGP solution suffers from ghost instability which makes the model 
unfavorable. Incorporation of a scalar field component on the DGP brane and treating 
the normal branch solutions brings a lots of new physics, some of which are studied 
in this paper. Previous studies in this field are restricted to either scalar fields on the 
self-accelerating branch or the simple case of minimally coupled scalar fields. In this 
paper we considered a scalar field component (both quintessence and phantom scalar 
fields), non-minimally coupled with induced gravity on the brane. We studied cosmo- 
logical dynamics of the normal branch solutions on the brane within a dynamical system 
approach. We translated dynamical equations into an autonomous dynamical system 
in each case. Then we obtained the critical points of the model in phase space of each 
model. The issue of stability of these solutions are studied with details. Also possibil- 
ity of having a stable attractor in de Sitter phase corresponding to current accelerated 
phase of universe expansion are studied in each case. We have also investigated the 
possibility to have a transition to the phantom phase of the equation of state parame- 
ter in each case. The classical stability of the solutions are treated also in a — w'^ 
phase-plane analysis in each step. We have shown in each step that there are several 
phases of accelerated expansion in each case, but only a limited critical points have 
stable, attractor solutions with de Sitter scale factor describing the current accelerated 
expansion on the brane. In summary, the main achievements of this study are as follows: 

• While the pure, normal DGP solution has not the potential to explain late-time 
cosmic acceleration, with a minimally coupled quintessence field in the normal DGP 
setup there is a stable de Sitter phase realizing the late-time cosmic speed up. Never- 
theless, as the pure DGP case, there is no possibility to cross the cosmological constant 
line by the effective equation of state parameter of the model. The classical stability 
domain of the model is restricted to those subspaces of the model parameter space that 
w'p < 3w v (l + Wp) with w v > —1 where w' v = ^f- and N = lna(t). 

• With a minimally coupled phantom field on the brane, it is possible to have an 
attractor, de Sitter solution realizing the late-time accelerated expansion in the nor- 
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mal DGP setup. Also, the effective equation of state parameter of the model crosses 
the phantom divide. Note that this crossing is impossible by the effective equation of 
state parameter of a minimally coupled quintessence field on the brane. Similar to the 
minimally quintessence field on the brane, the classical stability domain of the model 
is restricted to those subspaces of the model parameter space that w' < 3^(1 + w v ) 
with Wtp < — 1. 

• With a non-minimally coupled quintessence scalar field on the DGP brane, it is 
possible to realize a stable, de Sitter late-time accelerating phase in the normal branch 
of the model. It is possible also to cross the phantom divide by the effective equation of 
state parameter of the model. In this case, there are two different domains of stability 
in the w v — w' phase-plane of the model: a subspace with w'^ > 3w </ ,(l + w^) with 
Wtp < — 1 corresponding to an effective phantom phase and the other subspace with 
w' < 3u^(l + w v ) with w v > —1 corresponding to a quintessence phase. This feature 
shows that it is possible to have an effective phantom picture with a non-minimally 
quintessence field on the normal DGP setup. 

• For a non-minimally coupled phantom scalar field on the DGP brane, it is possible 
to realize a de Sitter late-time accelerating phase in the normal branch of the model. It is 
possible also to cross the phantom divide line by the effective equation of state parameter 
of the model in this case. Also, as for the case of non-minimally coupled quintessence 
field on the brane, there are two different domains of stability in the w v — w' phase- 
plane of the model: a subspace with w' > 310^,(1 + w^) with < — 1 corresponding 
to the phantom phase and the other subspace with w' < 3^(1 + w v ) with > — 1 
corresponding to an effective quintessence phase on the brane. Note that it is possible 
to have an effective quintessence picture with a non-minimally coupled phantom field 
on the normal DGP setup. 

Finally we stress that the observational status of the present DGP-inspired models 
can be treated in the same line as has been reported in Ref. [25] . 
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